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Abstract 

We consider a two-dimensional (d=2) kagome lattice gas model with attractive three-particle 
interactions around each triangular face of the kagome lattice. Exact solutions are obtained for 
multiparticle correlations along the liquid and vapor branches of the coexistence curve and at 
criticality. The correlation solutions are also determined along the continuation of the curvilinear 
diameter of the coexistence region into the disordered fluid region. The method generates a linear 
algebraic system of correlation identities with coefficients dependent only upon the interaction 
parameter. Using a priori knowledge of pertinent solutions for the density and elementary triplet 
correlation, one finds a closed and linearly independent set of correlation identities defined upon 
a spatially compact nine-site cluster of the kagome lattice. Resulting exact solution curves of the 
correlations are plotted and discussed as functions of the temperature, and are compared with 
corresponding results in a traditional kagome lattice gas having nearest-neighbor pair interactions. 
An example of application for the multiparticle correlations is demonstrated in cavitation theory. 
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I. INTRODUCTION 



Our fundamental understanding of anomalous thermodynamic and transport behaviors 
in a many-body cooperative system stems in large measure from concomitant knowledge 
and applications of its thermal equilibrium correlations. The familiar representations of 
macroscopic observables in terms of their underlying correlations, e.g., specific heat and 
magnetic susceptibility as energy and magnetization fluctuations, respectively, Kubo for- 
mulas, temperature-dependent Green's functions, fluctuation-dissipation theorems, and so 
forth, constitute many of the most instructive basic relationships and useful formulations 
in statistical physics and bestride virtually all areas of theoretical investigation in phase 
transitions, critical and multi-critical phenomena. Indeed, the impetus for the modern uni- 
fying interpretation of critical phenomena was the recognition of the essential role that the 
anomalously long-ranged spatial correlations played near a critical point resulting in scaling 
theories and the renormalization group approach towards problems in phase transitions and 
particle physics. Since spatial correlation functions are structured using thermal expecta- 
tion values of products of local variables, they clearly offer a more detailed description than 
thermodynamic for the order and symmetry in the system, and a precise presentation of the 
correlation solutions becomes highly desirable. 

Besides examining the asymptotically large distance behavior of correlation functions, it 
is also useful to obtain solutions for more spatially compact, short-distance type correlations 
These smaller scale correlations have varied applications, at criticality and otherwise. 
More particularly, Ising-type models have special appeal since, in select cases, their localized 
correlations [2] can be calculated exactly (any Hamiltonian having a finite density of finitely 
discrete commuting local variables can be cast as a standard or generalized Ising model). Ac- 
tually, the archetypical two-dimensional (d=2) Ising model [3] magnet with nearest-neighbor 
pair interactions and in zero magnetic field is the only realistic microscopic model of coop- 
erative phenomena for which many correlation solutions have been obtained exactly. Ising 
models are employed not only to represent certain kinds of highly anisotropic magnetic 
crystals but also, e.g., as lattice models for liquids, alloys, adsorbed monolayers, equilibrium 
polymerization, for biological and chemical systems, and in field theories of elementary par- 
ticles (lattice gauge theories describing the quark structure of hadrons). Applications for 
localized correlations in Ising-type models occur in the analysis of local equilibrium prop- 
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erties i„ the vicinity of isolated defects in the theory of both transport coefhcents ^ 

and thermodynamic response functions |7], in investigations of inelastic neutron scattering 
8] , percolation phenomena [ol, and in many other problems including topical connections 
between entanglements and spin correlations in quantum information theory [lOj], and a 
present example (Section VI) in cavitation theory 



Due to severe mathematical complexities, few analytically rigorous results are known for 
models having multiparticle interactions [12] ■ The lack of exact solutions for correlations in 
the models is an incentive for the present theoretical investigations. In the present paper, 
one obtains exact solutions for multiparticle correlations along the coexistence curve and 
at criticality of a kagome lattice gas with localized three-particle interactions. The phase 
diagrams for condensation of the lattice gas were determined previously [3] , specifically, the 
chemical potential and the density versus teinperature. The theory was based upon work by 
Wu [l^, Wu and Wu [ist, and Lin and Chen 16 1, who established that the partition function 
of a generalized (three-parameter) kagome Ising model having pair and triplet interactions 
and magnetic field is equivalent, aside from known pre-factors, to the partition function of 
a standard (two-parameter) honeycomb Ising model with pair interactions and field. Their 
theoretical developments incorporated a symmetric eight-vertex model on the honeycomb 
lattice in a mediating role. Later, the same result was obtained by Wu [17] using a direct 
mapping without a weak-graph transformation. 

Foreknowledge of the above phase diagrams (chemical potential and density vs tempera- 
ture) is vital supplementary information in the present theory, as is securing the solution for 
the elementary triplet correlation via appropriate logarithmic differentiation of the grand 
canonical partition function of the lattice gas. A linear algebraic system of correlation iden- 
tities is generated with coefficients dependent solely upon the interaction parameter of the 
model when evaluated along the coexistence curve. Using the known pertinent solutions 
for the density (no) and the elementary triplet correlation {nonin2) as a priori information, 
one succeeds in finding a closed and linearly independent set of correlation identities defined 
upon a spatially compact nine-site cluster of the kagome lattice. Resulting exact solution 
curves of the correlations are plotted and discussed as functions of the temperature, and 
are compared with the corresponding results in a traditional half-filled kagome lattice gas 
having nearest-neighbor pair interactions. To our knowledge, these are the first examples, 
certainly away from criticality, of exact solutions for multiparticle correlations in any pla- 
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nar lattice-statistical model with multiparticle interactions. Finally, the solutions for the 
multiparticle correlations are applied to cavitation theory in the condensation of the lattice 
gas. 

The paper is organized as follows. Section II presents the triplet-interaction kagome lat- 
tice gas model, and reviews its grand partition function and phase boundary curve. Section 
III reviews the liquid-vapor coexistence curve of the model, and calculates the relevant el- 
ementary triplet correlation. Section IV derives a basic generating equation for developing 
a linear algebraic system of correlation identities. Taking advantage of the supplemental 
information in Section III, Section V solves a set of the identities to secure solutions for 
multiparticle correlations along the coexistence curve and at criticality. Section VI provides 
an example of application for the multiparticle correlations in cavitation theory. Lastly, 
Section VII is a summary and discussion. 

II. PARTITION FUNCTION AND PHASE BOUNDARY CURVE OF A KAGOME 
LATTICE GAS WITH THREE PARTICLE INTERACTIONS 

Consider a lattice gas of TV atoms upon the kagome lattice (Figure 1) of jV sites with the 
(dimensionless) Hamiltonian 

-pn^K^, ^ riiUjUk, (2.1) 

<i,j,k> 

where f3 = l/ksT, with ks being the Boltzmann constant and T the absolute temperature, 
= j3e^ with 63 > being the strength parameter of the short-range attractive triplet in- 
teraction, the sum is taken over all elementary triangles, and the idempotent site occupation 
numbers are defined as 

{1 site I occupied, 
(2.2) 
site / empty. 

In (2.1), an infinitely-strong (hard core) repulsive pair potential has also been tacitly assumed 
for atoms on the same site, thereby preventing multiple occupancy of any site as reflected 
in the occupation numbers (2.2). 

In the usual context of the grand canonical ensemble, one introduces 

R = n-iiN (2.3) 
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where /x is the chemical potential with being the conjugate total number of particles 

N = Y,^i ■ (2-4) 

i 

Using (2.1), (2.3) and (2.4), the grand canonical partition function S(/i,A/', T) is given by 
S(/i,Ar,T) = ^e"^" = ^e^^^>""+^^^<'.^.*>"»"^"^ (2.5) 

where the summation symbol {rij} represents the set of total M occupation numbers. Aside 
from known pre-factors, the grand canonical partition function (2.5) on the kagome lattice 
can be transformed into the magnetic canonical partition function 

Z*{L*,K*) = ^e"'^^'- = J^e^*^'"'^^*^^'''^^'"', (2-6) 

{Mi} 

fii = ±1, i = I,-- - ,Af*, of a standard S = 1/2 Ising model ferromagnet upon the as- 
sociated honeycomb lattice (Figure 2) with (dimensionless) Hamiltonian —PTiJ^^ having a 
(dimensionless) external magnetic field L* and (dimensionless) nearest-neighbor pair inter- 
action parameter K* > 0, and where the summation symbol {/Xj} represents the set of total 
M*{= 2A/'/3) Ising variables. Specifically Q, 



<i,i,k> ■"■injTik 



el(^+M^\a* /2 cosh L*f K*)-i^* Z\L\K*), (2.7a) 



with the parameters 



[1 + (1 + e-/^^)-i/2]3/2 ' ^^■"') 
L* = ^ ln(e^3 -l)-^ ln(l + e"'^^), (2.7c) 
ir* = ^ln(l + e-^'^), (2.7d) 

and where the positivity K* > (ferromagnetic) is manifest in (2.7d). 

It is well known fisl] that a necessary and sufficient condition for the existence of a phase 
transition in the 5 = 1/2 ferromagnetic {K* > 0) honeycomb Ising model is the joint 
condition L* = and K* > K^, where the critical value K* = i ln(2 + ^3) = 0.65847 ■ • • . 
Imposing this joint condition of the associated Ising model upon the current calculations 
enables an exact solution to be found for the phase boundary curve of the triplet interaction 



kagome lattice gas model. In particular, the zero field condition (L* = 0) is realized by 
setting (2.7c) to zero which implies that the chemical potential is prescribed by the relation 

= (e^3 - 1)2/3 _ 1 at L* = 0. (2.8) 

Substituting (2.8) into (2.7d) gives 

= lln(e^« - 1) at L* = 0, (2.9) 
6 

relating the interaction parameters K*, K3 whenever the magnetic field parameter L* = 0. 
Using (2.9), the manifest positivity K* > from (2.7d) is therefore equivalent, at L* = 0, 
to In 2 < K3 < 00. In the next section, the restricted range ln2 < < 00 also assures 
that the solution found for the average particle number density p does not exceed unity 
(fully occupied lattice gas). For the remaining range < < In 2, the present model (2.1) 
doesn't admit phase transitions or critical behavior since the necessary zero-field {L* = 0) 
condition (2.8) cannot be satisfied by any real chemical potential fi. 

At criticality, the aforementioned literature value K* = ^ ln(2 + v^) is substituted into 
(2.9) yielding 

= ln[(2 + Vsf + 1] = 3.96992 ■ • • . (2.10) 

As comparison, for a traditional kagome lattice gas with attractive nearest-neighbor pair 
interactions, the corresponding critical value is known to be [sl, Q 

iTac = ln(3 + 2^3) = 1.86626- ■■ , (2.11) 

which is nearly 50% smaller than the critical value (2.10). The critical value K^c in (2.10) 
is used to locate the critical point in the relation (2.8), yielding the liquid-vapor phase 
boundary curve of the triplet interaction kagome lattice gas (Figure 3a). Explicitly, one 
directly obtains jisl 

/i/ea = -K^' ln[(e^^ - i f^ - 1], < K^JK, < 1, (2.12) 

with /i/e3 being a reduced chemical potential and K^clK^{= T/T^ a reduced tempera- 
ture where K^c{= es/keTc) = ln[(2 + ^3)3 + 1] = 3.96992 •■ • [(2.10)]. The curvilinear 
phase boundary begins at zero temperature with /i/es = —2/3 and ends (analytically) at 
a critical point whose coordinates are K^c/ = 1) A^/es = f^c/^s = —0.64469- ■ ■. At zero 
temperature, the phase boundary curve /x/es vs T/T^ has a zero slope in accordance with 



the Clausius-Clapeyron equation and third law of thermodynamics. Otherwise, its slope is 
positive which is more discernible at temperature closely below the critical temperature. As 
comparison, Figure 3b shows the corresponding phase diagram for the condensation of a 
conventional kagome lattice gas with attractive nearest-neighbor pair interactions. 



III. EXACT SOLUTIONS FOR (no), (noninz) ALONG COEXISTENCE CURVE OF 
KAGOME LATTICE GAS WITH THREE-PARTICLE INTERACTIONS 

The coexistence surface (and boundary) of the lattice-gas system is the prominent por- 
tion of its thermal equation of state surface. More particularly, the boundary (edge) of 
the coexistence surface contains the critical point and encloses the liquid-vapor coexistence 
region (heterogeneous mixture). The projective mapping of the boundary onto the chemical 
potential-temperature {fi — T) plane marks the previous liquid- vapor phase boundary curve 
(2.12), and a similar projective viewing in density-temperature (p — T) space indicates the 
liquid-vapor coexistence curve of the lattice gas. The latter phase diagram {p'[°y^ vs T /T^) 
: or the condensation of the triplet-interaction kagome lattice gas has already been obtained 
isl ]. A brief review of the results is now presented for later use. 

The solutions for the density p(= (no)) and the elementary triplet correlation (^0721^2) 
can be determined along the liquid- vapor coexistence curve and at criticality by logarithmic 
differentiations of the grand partition function S(/i,A/', T) [(2.7)] with respect to (3^ and 
i^3, respectively, and then letting L* for the full range of condensation temperatures 
< K^c/K^ < 1- Specifically, the exact solution for the average number density (no)(= p) 
is given, as L* — > 0, by 3] 

1 - ifl - (e^^ - + Wi)l.=o T < ^ < 1, 

1 - i[l - (e^^^ - l)-^/^][l + (/io/xi)L.=o], K % < = 5-72739 ■ ■ ■ , 

(3.1) 

where {poPi) l*=Oi {p)s are the nearest-neighbor pair correlation and spontaneous magneti- 
zation, respectively, of the previous 5* = 1/2 honeycomb Ising model ferromagnet. In the 
top of eq. (3.1), the — (+) algebraic sign corresponds to the hquid(vapor) branch of the 
coexistence curve. The restricted range ln2 < < 00 was required earlier to assure that 
the Ising interaction parameter K* > (ferromagnetic) and leads, using (2.10), to the fi- 
nite K^c/K^— range of the temperatures in the bottom of eq. (3.1). The terminating value 
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i^^3c/ln2 = 5.72739- • • also guarantees that p does not exceed pmax = 1 corresponding to 
the fully occupied ("close packed") lattice gas. The bottom expression of (3.1) is the con- 
tinuation of the curvilinear diameter of the coexistence region in the top expression of (3.1) 
beyond the coexistence surface. In (3.1), the exact solutions for {poPi) l*=o and {p)s are 
known to be Q, Q 

(/Xo/ii)L.=o = ^[coth2ir* + 7iri(K)], (3.2a) 

\{i-^^Y/\ o<§<i, 

{p)s = V -K - ^ ^3 

I 0, 1 < ^ < oo 

with Ki{k) being the complete elliptic integral of the first kind 

Ki{k)= (1- K^smH)-^/^de, (3.3a) 



and where 

= IQz^il + z^){l - zy^l - zY\ (3.3b) 

7 = (1 - z^){z^ -Az + - z^\{l - z)\ (3.3c) 

z = (3.3d) 

The expressions (3.1) can be written purely in the natural notation of the lattice 
gas by substituting the interaction parameter relation (2.9) into (3.2) and (3.3). Then, 
substituting the resulting forms into the top expression of (3.1) yields the sought exact 
solution for the liquid-vapor coexistence curve which is plotted in Figure 4a. Additionally, 
the composition of the density p at points within the two-phase coexistence region in Figure 



4a can be determined by the application of the thermodynamic lever rule 2l|]. The phase 
diagram of Figure 4a exhibits an asymmetric rounded shape with a curvilinear diameter, 
contrasting the familiar symmetric rounded shape and constant rectilinear diameter for a 
conventional kagome lattice gas with attractive nearest-neighbor pair interactions (see Figure 
4b). In the global comparison (i.e., "over- laying") of Figure 4a and Figure 4b, the values 
of the liquid-vapor branches (and curvilinear diameter) in Figure 4a exceed the values of 
the liquid-vapor branches (and rectilinear diameter) in Figure 4b at all corresponding finite 
reduced temperatures. 

The curvilinear diameter of the coexistence region in Figure 4a is the solid curve which 
begins at zero temperature with p = 1/2, exhibits a positive slope which is more pronounced 
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at temperatures closely below the critical temperature, and ends at the critical point (solid 
circle) whose coordinates are Ksc/Xs — \^ p — p^ — 0.58931 • • • (see (3.7a)). Using the bot- 
tom expression of (3.1), the curvilinear diameter is continued beyond the coexistence region 
as the solid sigmoidal curve, eventually ending at the point (solid square) with coordinates 
K^JK^ = K^J In 2 = 5.72739 ■ ■ • , p = p„„. = 1. 

To obtain the exact solution for the triplet correlation (nonin2) along the coexistence 
curve, one first computes the logarithmic derivative of S(/x, jV, T) with respect to the triplet- 
interaction parameter K^. Specifically, (2.7a) gives 

lnS(//,A/',r) = + 13 p)M* a* -\u cosh L* -\u2) 

- ^M*hicoshK* ^hiZ*{L\K*), (3.4) 

where the total number of kagome lattice sites M — \M* with M* being, as stated previously, 
the total number of lattice sites of the associated honeycomb lattice (see Figure 2). Also, 
the total number of elementary triangles A/a of the kagome lattice equals A''*, again seen in 
Figure 2. Using (2.6), (2.7) and (3.4), one obtains 

ainS(At,AA, T) 



"^^"«*-(tanhr)-^^^ 



hA^ , ...^ , 

3*r*/. ^\.*^^K* din Z* dL* ^ din Z*dK* ,^ , 

-AA*(tanhX*)^ + + (3.5a) 



where 



dlnZ* 
dL* 

dL* 1 e^' 



= (3.5b) 

dK, 2 e^3 _ 1 ' (3.5c) 

OK* , „ 

= 0, (3.5d) 



dK- 
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aina* 1 



dK^ 8 (e-/^^ + 1)3/2 + e^3 _ 1 

Hence, as L* — > 0, (3.5) yields 

1 

(nonin2) = - ^^^^ _ ^ (1 ± ky)s), ln2 < < oo 

life(l±(/^>5), 0<^<1, 
li^' K ^ = 5.72739- 



(3.6) 



having replaced (/ij) by ±(/i)s in (3.5b) due to 'spontaneous symmetry breaking', hav- 
ing substituted the zero-field (L* = 0) condition (2.8) into (3.5e), used the identification 
A/a = A/"*, and having partitioned, as earlier, the (dimensionless) inverse temperature range 
ln2 < < oo into (reduced) temperature intervals below and above criticality in (3.6). 
Using the known result (3.2b) for the spontaneous magnetization (/i)5, and substituting the 
interaction parameter relation (2.9) into (3.2b) and (3.3b,d), the exact solution in the up- 
per expression of (3.6) for the elementary triplet correlation (^0^1722) is obtained along the 
coexistence curve and at criticality as a function of the (reduced) temperature K^d and 
is plotted in Figure 5a. Also shown in Figure 5a, the lower expression of (3.6) is the exact 
solution for (rio'^i^2) along the continuation of the curvilinear diameter (solid curve) beyond 
the coexistence surface, where the trajectory terminates at the point (solid square) with co- 
ordinates K-ic/K^ = 5.72739- ■■, {nQnin2) = 1. The corresponding results for (^0^1712) in 
a kagome lattice gas with attractive nearest-neighbor pair interactions are shown in Figure 
5b where, as a function of temperature, the curvilinear diameter and its extension into the 
disordered fiuid region exhibit a sigmoidal shape and monotonically decreasing behavior 
contrasting the results shown in Figure 5a for the three-particle interactions. 



At criticality, the honeycomb Ising ferromagnet has values {fi)s = and [23] {poPi)c = 
A\/3/9. These critical values along with the critical value (2.10) e^^" = (2 + v^)^ + 1 are 
substituted into (3.1) and (3.6) yielding, respectively, 

Pc = ^(^-V3) = 0.58931- ■• , (3.7a) 

{non,n2)c = ^( ^^^j^^ ) = 0.50961 • ■ ■ , (3.7b) 

as critical values in the triplet-interaction kagome lattice gas. For a traditional kagome 
lattice gas with attractive nearest-neighbor pair interactions, the corresponding values are 

Pc = ^ = 0.50000- ■■ , (3.8a) 

(nonin2)c = — (3 + 2^3) = 0.40400 ■ ■ ■ . (3.8b) 
16 

As seen, the one-parameter (es) model Hamiltonian (2.1) exhibits a single critical point on 
its thermal equation of state surface ^ = fi{p, T) (in suitable units), with critical coordinates 
^Je^ = -0.64469 ■ ■ ■ , pc = 0.58931 ■ ■ ■ , kBTje^ = K^^ = 0.25189 ■ ■ ■ . The coexistence 
surface, its boundary coexistence curve and the critical point comprise the prominent portion 
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of the above equation of state surface. The hquid and vapor branches of the coexistence 
curve identify at the critical point, and all lattice-gas thermal averages in the current studies 
are evaluated along both branches and at criticality. 

Some comments are warranted on the nature of the mathematical singularities in the 
phase diagram of Figure 4a. The critical behaviors of the coexistence curve and the curvi- 
linear diameter are underlaid by the known critical behaviors of the honeycomb Ising thermal 
averages (/i)s and (/io/^i) l*=o in (3.1). Letting pi, denote the particle number density along 
the liquid and vapor branches, respectively, of the coexistence curve, the top expression in 
(3.1) immediately reveals that the ordering parameter pi — p^ (length of vertical "tie-line" 
spanning the coexistence region) realizes its critical behavior purely from {p)s, whereas the 
curvilinear diameter |(p; + p^) (arithmetic mean of pi and p„) realizes its critical behavior 
solely from {poPi) l'=o- Hence, the relation (3.2b) for the Ising spontaneous magnetization 
{p)s leads to the vanishing of the order parameter pi — p^ at the critical point with the Ising- 
type critical exponent 1/8 (algebraic branch point singularity). Similarly, relation (3.2a) for 
the Ising nearest-neighbor pair correlation (popi) l*=o leads to the result that the curvilin- 
ear diameter ^{pi + p^) and its analytic continuation into the disordered fluid region (lower 



expression of (3.1)) possess a weak Ising energy-type singularity 23| elne at the critical 
point, where e > is a small fractional deviation of the temperature from its critical value. 
In contrast, note that the constant rectilinear diameter in a conventional d = 2 lattice gas 
(Figure 4b) is analytic at criticality. One also recognizes that the curvilinear diameter of 
the solution curve for the elementary triplet correlation in (3-6) and Figure 5a is 

analytic at the critical point. In the above mathematical arguments, one uses (2.9), (2.10) 
and the previous literature value K* = | ln(2 + -\/3) to establish an exact scaling relation 

between the smallness parameters e* and e. The scaling relation (3.9) affords a direct proof 
that the singularities in the associated Ising model lead to the same nature of singularities 
in the lattice gas phase diagram in (3.1) and Figure 4a, viz., e*^/^ — >■ e^^^, e*lne* elne, 
where, neglecting second-order small quantities, e* = {K* — K*)/K*, e = {K^ — K2,c)/K2,c- 
In the present section, we emphasized that the solutions for (no) and {nQnin2) were 
determined along the coexistence curve via logarithmic differentiations of the grand partition 
function E{p,J\f,T) with respect to Pfi and K^, respectively, and then letting L* 0. To 
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secure solutions for additional correlations, a different method is needed. In the subsequent 
sections, one develops and solves, as L* 0, a linear algebraic system of correlation identities 
whose coefficients depend solely upon the interaction parameter K^. The a priori knowledge 
of (jiq) and {nQnin2) will be instrumental in the quest for closure and linear independence 
within the linear algebraic system of identities. 



IV. BASIC GENERATING EQUATION FOR CORRELATION IDENTITIES 

The class of correlation identities currently considered is a set of linear algebraic equations 



24| . To develop such 



with coefficients dependent only upon the interaction parameter K3 
identities systematically, one proceeds to derive their basic generating equation. 

Let [g] be any function of the lattice-gas variables ni,n2,--- ,nj^_i {excluding uq, the 
origin site variable in Figure 1). Letting H' = Ti' — j^N', Yl'{n } denote a restricted energy 
form and summation operation, respectively, which exclude no, one can construct the grand 
canonical thermal average {nQ[g]) as 

{ui} no 



{Ui} 



j^^Qno[f3fj.+K3{nin2+n-in4)] 



EpUo [/311+K3 (nin2 +n3n4)] 



(4.1) 



yielding 



^f3fi+Ki(n-in2+n3n4) 

K[^]) = ( ^^gff^+^3(^,^,+^3^,) [^?]), no<^[g], (4.2) 



having lastly used the standard definition of grand canonical thermal average initiating (4.1). 
In deriving (4.1), the "split, rearrange, then reconstitute" procedures are justified since all 
lattice-gas variables commute. To further develop (4.2), one writes 

p/3/i+ii'3(nin2+n4n4) 

A + _B(nin2 + 713^4) + Cnin2n3n4, (4.3) 



where the expansion as a finite algebraic series in the lattice-gas product variables ni?7,2, ^377,4 
reflects their idempotent nature [(npn^)^ = npUg]. The coefficients A, B, C are determined 
by considering the following realizations of the product variables nin2, n^n^^ in (4.3): 
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(i) nin2 — rizUi — 0, yielding 



^ = IT = 7' (4-4a) 

1 + ef^f" e-f^^" + 1 ^ ^ 



(ii) 711^2 = 0, ^3714 = 1, yielding 



+ ^ = T 

(iii) nin2 = n^n^ = 1, yielding 



A + B = -r—-r- = — 3 — , (4.4b) 



/3/i+2j£r3 ]^ 

The three linear algebraic inhomogeneous equations (4.4 a-c) in the three unknowns A, B, 
C are patently linearly-independent and directly give the solutions 

A = / ^ , (4.5a) 

B = . \ -^——^ (4.5b) 

1 2 1 

For adoption along the coexistence curve, one evaluates the coefficients (4.5) at L* — 0. 
Substituting the (2.8) expression into (4.5), one obtains, as L* —>■ 0, 

A = (e^3 - l)-'/^ (4.6a) 

B = e^n(e^^ - 1)2/3 + - 1]-^ - (e^^ - l)-^/^, (4.6b) 
C = e2^3[(ei^3_ 1)2/3 ^g2/f3_;^]-i 

_ 2e^^ [(e^^ - 1)2/3 + - 1]-^ + (e-^^ - l)-^/^. (4.6c) 

Returning to (4.2) and substituting (4.3), the basic generating equation for multiparticle 
correlation identities is thus given by 

{no[g]) = Mig]) + B{{nin2 + n3n4)[g]) + (7(^1712713^4 [gr]), no ^ [5-], (4.7) 

where the coefficients A, B, C are given by (4.6) as L* 0. The linear algebraic system of 
correlation identities generated by (4.7) will be employed in the next section to determine 
the solutions for various multiparticle correlations along the coexistence curve of the triplet- 
interaction kagome lattice gas. 
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V. EXACT SOLUTIONS FOR CORRELATIONS ALONG THE COEXISTENCE 
CURVE OF KAGOME LATTICE GAS WITH THREE-PARTICLE INTERAC- 
TIONS. 

Solutions will now be determined for correlations along the coexistence curve of the 
present lattice gas model. More particularly, exact solutions arc found for the nearest- 
neighbor pair and various multiparticle correlations defined upon the select nine-site cluster 
shown in Figure 6. In section III, pertinent exact solutions were obtained for (0) and (012), 
where, for notational simplicity, only the numeric site labels in Figure 6 are written inside 
the thermal average symbols. As seen shortly, the a priori knowledge of (0) and (012) is 
vital supplementary information within the system of correlation identities. 

First consider the five-site ("bow-tie") cluster 0, 1,2,3,4 (see Figure 6) and the five gen- 
erators (0), (01), (012), (0123), (01234). Employing the basic generating equation (4.7), one 
directly obtains the following correlation identities for (0) and (012), respectively, 

(0) = ^ + 5(12 + 34) + C(1234), (5.1a) 
(012) = A(12) + 5(12 + 1234) + C(1234), (5.1b) 

having used the idempotent property n| = rij. The identity (5.1a) is initially the only 
inhomogeneous equation in the infinite system of identities. However, since (0) and (012) 
are already known along the coexistence curve, one rearranges (5.1) into the standard form 
of two linear algebraic inhomogeneous equations in the two unknowns (01) and (1234): 

25(01) + C(1234) = (0) - A, (5.2a) 

(A + 5)(01) + (5 + C)(1234) = (012), (5.2b) 

where the symmetry of the kagome lattice has been recognized in the equating of 
geometrically-equivalent pair correlations in (5.2). Since the LHS coefficient matrix in (5.2) 
has a non-vanishing determinant, equations (5.2) arc linearly independent and hence deter- 
mine the solutions for the nearest-neighbor pair (01) and quartet (1234) correlations along 
the coexistence curve. 

Again employing the basic generating equation (4.7), the three remaining generators (01), 



14 



(0123) and (01234) yield, respectively, the identities 

(01) = A{1) + B{12 + 134) + C(1234), (5.3a) 

(0123) = (>1 + S) (123) + (S + C) (1234), (5.3b) 

(01234) = {A + 2B + C){12U). (5.3c) 



The relevant exact solutions for the thermal averages (0), (01), (1234) appearing in (5.3) have 
already been obtained. Thus, using lattice symmetry ((1) = (0), (12) = (01), (134) = (123)), 
the identity(5.3a) directly determines the triplet correlation (123), identity (5.3b) in turn 
directly establishes the quartet generator (0123) since the RHS is known, and identity (5.3c) 
similarly secures the quintet generator (01234) since the RHS is known. To review, exact 
solutions are now known along the coexistence curve for the density (0) and six correlations 
defined upon the five-site ("bow tie") cluster 0, 1, 2, 3, 4 in Figure 6, specifically, 

(0), (01), (012), (123), (0123), (1234), (01234). (5.4) 

One next considers the seven-site cluster 0, 1, 2, 3, 4, 5, 6 (see Figure 6) and the six generators 
< 03456 >, < 0456 >, < 01256 >, (012356), (012456), (0123456). Using the basic generating 
equation (4.7), the latter generators yield, respectively, the identities 



(03456) 


= {A + B) (3456) + {B + C) (123456) , 


(5.5a) 


(0456) 


= A(456) + S(3456 + 12456) + C(123456), 


(5.5b) 


(01256) 


= (A + fi)(1256) + (5 + C)(123456), 


(5.5c) 


(012356) 


= (A + 5)(12356) + (B + C) (123456), 


(5.5d) 


(012456) 


= {A + B) (12456) + (S + C) (123456) , 


(5.5e) 


(0123456) 


= (A + 25 + C) (123456). 


(5.5f) 



Using the lattice symmetry, some correlations in identities (5.5) on the seven-site cluster are 
in registry with known correlations (5.4) on the earlier five-site cluster. Specifically, 

(03456) = (01234), (3456) = (0456) = (0123), (456) = (012). (5.6) 

The registry listings (5.6) are now used as a priori information in identities (5.5) thereby 
reducing the number of unknown correlations. Hence, in identity (5.5a), all correlations 
are known except the sextet correlation (123456)(= (012356)), so the latter is determined. 
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Similarly, in identity (5.5b), all correlations are now known except the quintet correlation 
(12456), so the latter is obtained. In identity (5.5c), all correlations are known save the 
quartet correlation (1256), so the latter is secured. Continuing the cascade reasoning, the 
only unknown correlation in (5.5d) is the quintet correlation (12356) so the latter is found. In 
identity (5.5e), all correlations are known save the sextet generator (012456), so the latter 
is determined. Lastly, in identity (5.5f), the RHS correlation is known so the LHS septet 
generator (0123456) is obtained. In summary, exact solutions have been found along the 
coexistence curve for the additional six correlations defined upon the seven-site cluster 
0,1, 2, 3, 4, 5, 6 in Figure 6: 

(1256), (12356), (12456), (123456), (012456), (0123456). (5.7) 

One proceeds to consider the nine-site cluster 0, 1, 2, 3, 4, 5, 6, 7, 8 (see Figure 6) and the six 
generators (0345678), (012345678), (05678), (045678), (0125678), (01245678). Employing 
again the basic generating equation (4.7), the above generators yield, respectively, 

(0345678) = (A + S) (345678) + (S + C) (12345678), (5.8a) 

(012345678) = (A + 25 + C)( 12345678), (5.8b) 

(05678) = ^(5678) + 5(125678 + 345678) + C(12345678), (5.8c) 

(045678) = A(45678) + 5(1245678 + 345678) + C(12345678), (5.8d) 

(0125678) = (A + 5)(125678) + (5 + C)(12345678), (5.8e) 

(01245678) = (^ + 5)(1245678) + (5 + C) (12345678). (5.8f) 

Using lattice symmetry, various correlations in (5.8) on the nine-site cluster are in registry 
with known correlations (5.7) on the previous seven-site cluster. Specifically, 

(0345678) = (0123456), (345678) = (012456), (05678) = (12356), 

(5678) = (1256), (045678) = (123456), (45678) = (12456). (5.9) 

The registry listings (5.9) are now used as a priori information in (5.8) thus reducing the 
number of unknown correlations and enabling one to elicit further correlation solutions from 
the system of identities. In identity (5.8a), all correlations are known except the octet 
correlation (12345678) so the latter is determined. Hence, the RHS of identity (5.8b) is 
estabhshed thereby yielding the LHS nonuplet generator (012345678). In identity (5.8c), 
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all correlations are known save the sextet correlation (125678) so the latter is obtained. In 
identity (5.8d), all correlations are known save the septet correlation (1245678) thus the 
latter is secured. In identity (5.8e), both RHS correlations are known so the LHS septet 
generator (0125678) is found. Lastly, in identity (5.8f), both RHS correlations are known, 
hence the LHS octet generator (01245678) is determined. In summary, exact solutions have 
been obtained along the coexistence curve for the additional six correlations defined upon 
the nine-site cluster 0, 1, 2, 3, 4, 5, 6, 7, 8 in Figure 6: 

(125678), (1245678), (01245678), (0125678), (12345678), (012345678). (5.10) 

To review, one considered correlations defined upon the nine-site cluster (Figure 6) of the 
kagome lattice comprising a central triangle with its three corner-sharing triangles. The col- 
lective contents of (5.4), (5.7) and (5.10) reveal that exact solutions have been found for the 
density (0), nearest-neighbor pair correlation (01), and seventeen multiparticle correlations, 
all along the coexistence curve and at criticality, of the triplet-interaction kagome lattice 
gas. Also, the solution for each thermal average was determined along a finite-length critical 
isofield {L* = 0) trajectory in the disordered fluid region of p — T space. In this context, 
illustrative correlations (01), (1234), (01234), (012345678) are each plotted as functions of 
the reduced temperature x — K^c/ K3{— T/T^ in Figures 7a,8,9,10, respectively. Along 
the coexistence curve (Figure 4a), each correlation solution exhibits an asymmetric rounded 
shape with an infinite slope at the critical point (solid circle). The degree of asymmetry is 
associated with the slope of the curvilinear diameter in the graph, which is more pronounced 
at temperatures closely below the critical temperature (x = 1). One observes that the slope 
is positive for the elementary triplet and nearest-neighbor pair correlations in Figures 5a and 
7a, respectively, and negative for the quartet, quintet and nonuplet correlations in Figures 
8, 9 and 10, respectively. A positive (negative) slope of the curvilinear diameter implies 
that the correlation solution curve along the upper liquid branch falls slower (faster) than 
the corresponding rise of the solution curve along the lower vapor branch. In each case of 
negative slope of the curvilinear diameter (Figures 8, 9, 10), the correlation solution along 
the continuation of the curvilinear diameter into the disordered region shows a rounded 
minimum before ascending to the common maximum value (solid square) associated with 
the fully occupied lattice. 
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VI. EXAMPLE OF APPLICATION 



A simple example employing multiparticle correlations occurs in cavitation theory 
One seeks the probability of an elementary cavity or void in the fluid system, where an 
elementary cavity is construed to be an empty site with all its nearest-neighbor sites occu- 
pied. In a kagome lattice gas, the joint probability po,!, 1,1,1 ^^^^ origin site is empty 
and its four nearest-neighbor sites 1,2,3,4 are occupied equals the thermal average value 
((1 — nQ)nin2n^ni) . This probabilistic interpretation of the thermal average can be directly 
extracted from the standard deflnition of average value in probability theory, as will now be 
shown. 

Let /(no, rii..., n/s) be a function of the lattice-gas discrete random variables ng, ni, ...,77,4 
(see "bow-tie" cluster in Figure 1). Then, the average value (/(rio, ni, ...,77,4)) is deflned by 

Pno,ni,...,n4 

/(rao,ni, ...,774), (6.1) 

{no,ni,...,n4} 

where Pno,ni,...,n4 is the joint probability of a specific realization of the variables no, ni, ...,714 
and the summation symbol {nQ,ni, ...,■^,4} represents the set of all possible 2^ = 32 realiza- 
tions of the flve lattice gas variables. 
Considering a select function 

f{no, ni, 124) = (1 - no)nira2n3n4, (6.2) 

the deflnition (16.11) immediately "fllters out" the desired result 

((1 - no)nin2nsni) = Po,i, 1,1,1, (6.3) 

since the value of /(0,1,1,1,1) = 1 and all other values of the function (16. 2p vanish in the 
weighted summation (16.11) . This completes the proof. 

It is instructive to mention that taking thermal averages of select operator functions is a 
direct and general method for exactly structuring the elements of reduced statistical density 
matrices [25j in equilibrium statistical mechanics. Each diagonal element of a reduced statis- 
tical density matrix is a thermal average value having a familiar probabilistic interpretation 
similar to (16.31) . 

Expanding (16.31) . the probability po,i,i,i,i(= p) of an elementary cavity becomes the dif- 
ference 

P = {nin2n2,ni) - {nQnin2n2,n^) . (6.4) 
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The representation (16 ■4p is valid for arbitrary lattice gas interactions (pair, triplet, ...) upon 
any four-coordinated regular lattice, i.e., the d = 2 square and kagome lattices and the 
d = 3 diamond lattice. In Section V, exact solutions for both RHS correlations in (16. 4p 
were obtained along the coexistence curve and at the critical point of the triplet-interaction 
kagome lattice gas. Consequently, the probability p in (16.41) is likewise determined along the 
coexistence curve and at criticality, as well as along the finite-length critical isofield {L* = 0) 
curve in the disordered fluid region of p — T space. In this setting, the RHS difference of 
correlations in (6.4) is seen as the graphical subtraction of Figure 9 from Figure 8. 

The results, viz., p vs x, are exhibited in Figure 11 and Figure 12, the latter being an 
enlargement of the former in the range of (reduced) condensation temperatures < x < 1. 
At absolute zero temperature, the coexistence curve p vs a; in Figure 4a is unity (zero) 
for the liquid (vapor) branch. In either case (fully occupied or completely empty lattice), 
an elementary cavity cannot exist, so p vanishes at zero temperature in Figure 11 (or 12). 
For increasing temperatures in Figure 12, the probability p of an elementary cavity along 
the liquid branch increasingly exceeds that along the vapor branch, except at temperatures 
closely below the critical temperature {x = 1) where the two branches become identical. Such 
asymmetric behavior is physically anticipated since decreasing the particle number density 
pi along the hquid branch, beginning from a fully-occupied ground state liquid phase, creates 
increasingly more elementary cavities than increasing py along the vapor branch, beginning 
from a completely empty ground state vapor phase. Furthermore, the curve p vs x in the 
disordered region {x > 1) of Figure 11 exhibits a rounded maximum Pmax = 0.0409 ■ ■ • 
at X = X = 3.6 ■ ■ ■ , and eventually vanishes at a; = 5.72739 ■ ■ ■ corresponding to p = 
Pmax = 1 (fully occupied lattice) in Figure 4a. The latter node of p along the critical isofield 
[L* = 0) trajectory is deemed to be a multiparticle interaction effect. In Figure 11, the 
above maximum value Pmax is larger by an order of magnitude than the lower temperature 
(x = 0.98 ■ ■ ■ ) maximum value along the liquid branch of the coexistence curve. A synoptic 
view of Figures 4a, 11 and 12 enables one to numerically eliminate their common temperature 
variable x and obtain the graph p vs p. In particular, the probability value Pmax = 0.0409 ■ ■ • 
corresponds to a density value p = p = 0.88 ■ ■ ■ , and p monotonically increases (decreases) 
for the density range pc < p < p{p < p < 1), where pc = 0.58931 • ■ ■ [(3.7a)]. 

Reflecting upon the traditional case of attractive nearest-neighbor pair interactions in a 
kagome lattice gas, the critical isofield (L* = 0) solution for p in the disordered region of 
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Figure 4b is a constant at all temperatures, i.e., 1 < x < oo, p = | {half-filled lattice). The 
half-filled condition infers that the probability p of an elementary cavity is non-vanishing 
along the critical density (pc = |) line at all temperatures 1 < x < oo in Figure 4b. For 
instance, as the (reduced) temperature x — oo, p — > (i)^ — (i)^ = M ~ 0.03125, having 
used (16.41) with the realization therein that the average of a product equals the product of 
the averages at infinite temperature. The absence of a node in the solution curve p vs x for 
pair interactions at all temperatures 1 < x < oo lends credence to the earlier contention that 
the finite-temperature node in Figure 11 for triplet interactions is a multiparticle interaction 
effect. 



VII. SUMMARY AND DISCUSSION 

Exact results in physics are valuable for a variety of reasons. Endeavoring to retain 
only the most essential ingredients of a physical problem, exact solutions of simple model 
systems often provide definite guidance and insights on more realistic and invariably more 
mathematically complex systems. Exact results in tractable models of seemingly different 
physical systems may alert researchers to significant common features of these systems and 
actually emphasize concepts of universality. In addition to their own aesthetic appeal, exact 
results can, of course, serve as standards against which both approximation methods and 
approximate results may be appraised. Also, the underlying mathematical structures of 
exactly soluble models in statistical physics are rich in content and have led to important 
developments in mathematics. 

Phase diagrams for the condensation of a two-dimensional (d=2) kagome lattice gas with 
purely three-particle interactions were determined previously 13|. Specifically, the liquid- 
vapor phase boundary (chemical potential vs temperature) and the companion liquid- vapor 
coexistence curve (density vs temperature) were obtained. These exact phase diagrams 
were briefly reviewed in the present paper to make it self-contained and for later key use 
in securing solutions for multiparticle correlations along the coexistence curve and at the 
critical point of the triplet-interaction kagome lattice gas. Pertinent exact solutions for 
the average particle-number density (uq) and the elementary triplet correlation {nonin2) of 
the lattice gas were obtained via direct logarithmic differentiations of the grand canonical 
partition function with respect to the (dimensionless) variables /5p and K^, respectively, and 
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then taking an appropriate vanishing field hmit (L* 0). However, to estabhsh solutions 
for additional correlations upon the coexistence curve, a different method was necessary. 
The method required supplemental use of the known solutions for (ng) and (^0^177,2), and is 
briefly discussed below. 

A linear algebraic system of correlation identities was generated having coefficients depen- 
dent solely upon the interaction parameter when evaluated along the coexistence curve. 
Using the previously determined solutions for (uq) and {nQnin2) as a priori information, one 
succeeded in finding a closed and linearly independent set of correlation identities defined 
upon a spatially compact nine-site cluster (Figure 6) of the kagome lattice. Employing sim- 
ple algebraic techniques, exact solutions are thereby now known for the density (no), the 
nearest-neighbor pair correlation {nQiii) and seventeen multiparticle correlations, all along 
the coexistence curve and at criticality of the triplet interaction kagome lattice gas. The 
correlation solutions were also secured along the continuation of the curvilinear diameter of 
the coexistence region into the disordered fiuid region of p — T space. To our knowledge, 
these are the first examples, certainly away from criticality, of exact solutions for multiparti- 
cle correlations in any planar lattice-statistical model with multiparticle interactions. Along 
the coexistence curve (Figure 4a), the graphs of various correlation solutions were plotted 
(Figures 5a, 7a, 8, 9, 10) as functions of the temperature, and each solution curve exhib- 
ited an asymmetric rounded shape with an infinite slope at the critical point (solid circle). 
The degree of asymmetry was associated with the slope of the curvilinear diameter in the 
graph, and a negative slope (Figures 8, 9, 10) was a harbinger for the correlation solution 
to have a rounded minimum along the continuation of the curvilinear diameter into the 
disordered region. The exact solution curves were compared with the corresponding results 
in a conventional half-filled kagome lattice gas having nearest-neighbor pair interactions, 
and the solutions for the multiparticle correlations were applied to cavitation theory in the 
condensation of the lattice gas. 

The class of correlation identities currently considered has appeared in the literature 24 1 
for almost five decades, particularly in the context of planar Ising models with nearest- 
neighbor pair interactions. In the present paper, however, the correlation identities were 
developed within the framework of the grand canonical ensemble for a kagome lattice gas 
having localized three-particle interactions. As reviewed above, solutions were determined 
for various localized correlations of the lattice gas. One further reasons that the lengths of 



21 



the vertical "tie-lines" connecting the upper and lower branches of the solution curves vanish 
with an Ising-type critical exponent 1/8, and the curvilinear diameter of each solution curve 
is either analytic at the critical point or possesses a weak Ising energy-type singularity e In e. 
Supporting arguments for these asserted behaviors involve the following particulars. The 
critical behaviors of the lattice gas model are embedded in the system of identities at the 
outset, entering the initial identities (5.2a,b) via substitution of known solutions for (no), 
(^0/11^2) into the inhomogeneous terms. The remaining identities (5.3), (5.5) and (5.8) are 
homogeneous relations among correlations. Solutions for seventeen correlations are then 
directly determined by elementary linear algebraic procedures and cascade orderings within 
the system of identities. Aside from well-behaved additive terms containing solely coefficients 
(4.6), all solutions can eventually be represented as linear combinations of (no) and (nonin2). 
Hence, one concludes that the familiar Ising-type critical singularities e^/®, elne inherent in 
(no) and (nonin2) similarly characterize the singular behaviors in the subsequent seventeen 
correlation solutions. 

Extensions in these results for the number of correlation solutions and the shapes and 
sizes of the n-site clusters could be aided by computer search codes in linear algebra theory. 
In these searches, one anticipates linear independence in a system of identities to be a more 
elusive algebraic property than closure. Finally, if one assumes a ''clustering property" 
(asymptotic separability) for the correlations in the model, one can then employ the present 
known solutions for correlations on smaller length scales to gain information on chosen 
asymptotically long-distance-type correlations. A simple example reifies these concepts. Let 
Qi, 92; ri, r2, rs; si, S2, S3, S4, S5 be the lattice sites of a nearest-neighbor bond, an elementary 
triangle and a "bow-tie" configuration, respectively, where the q—, r— , s— clusters of sites 
are mutually separated by asymptotically large distances. Then, the "clustering property" 
yields the decuplet correlation 

)- (01) (012) (01234), (7.1) 

having lastly used lattice symmetry and earlier notations. The result (7.1) shows that the 
chosen decuplet correlation asymptotically equals the graphical multiplication of the solution 
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curves in Figures 5a, 7a and 9. 
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FIG. 1: The kagome lattice is a two-dimensional periodic array of equilateral triangles and regular 
hexagons. The lattice is regular (all sites equivalent, all bonds equivalent) and has coordination 
number 4. The origin site and its four nearest neighboring sites are specifically enumerated. 
Whenever three atoms of the lattice gas simultaneously occupy the vertices of an elementary 
triangle (say sites 0, 1,2), these atoms experience a triplet interaction with strength parameter 63. 




FIG. 2: A honeycomb lattice (dashed bonds) may be associated with the kagome lattice (solid 
bonds). A honeycomb lattice (dashed bonds) is a two-dimensional periodic array of regular 

hexagons, is regular (all sites equivalent, all bonds equivalent) and has coordination number 3. 
The origin site and its three nearest-neighboring sites are specifically enumerated upon the associ- 
ated honeycomb lattice. 
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FIG. 3: (a) The liquid-vapor phase boundary of the triplet-interaction kagome lattice gas, with 
///es being a reduced chemical potential and x = K^d K^{= T/Tc) a reduced temperature, where 
K^c{= ^a/ksTc) = ln[(2 + a/S)^ + 1] = 3.96992- • • . The curvilinear phase boundary begins at 
zero temperature with /u/ea = —2/3, and ends (analytically) at the critical point (solid circle) 
whose coordinates are x = 1, n/es = njes = - ln(4A/3 + 6)/ ln(15V3 + 27) = -0.64469 • • • . The 
continuation of the phase boundary is the dashed curve which eventually diverges logarithmically 
at the reduced temperature x = i^3c/ln2 = 5.72739 ■■■ . (b) The liquid- vapor phase boundary 
of a conventional kagome lattice gas having nearest-neighbor pair interactions, with /x/e2 being a 
reduced chemical potential where £2 > is the strength parameter of the attractive pair interaction, 
and X = K2c/K2{= T/Tc) being a reduced temperature, where K2c = ln[(3 -|- 2\/3)] = 1.86626 • • • . 
The phase boundary is rectilinear, specially, the constant valued expression /x/e2 = —2, < a; < 1 
(solid circle locates the critical point). The continuation of the phase boundary is the horizontal 
dashed line which continues indefinitely. 
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FIG. 4: (a) The liquid-vapor coexistence curve of the triplet-interaction kagome lattice gas, with p 
being the particle number density and x = K^c/K-i{= T/Tc) a reduced temperature, where K^d^ 
es/ksTc) = ln[(2 -|- \/3)^ + 1] = 3.96992 • • • . The curvilinear diameter of the asymmetric rounded 
coexistence region is the solid curve which begins at zero temperature with p = 1/2, and ends at the 
critical point (solid circle) whose coordinates are x = I, p = pc = |(7/2 — \/3) = 0.58931 • • • . The 
continuation of the curvilinear diameter into the disordered fluid region is the continuing sigmoidal 
curve which is monotonically increasing and concave downward, eventually ending at the point 
(solid square) with coordinates x = i('3c/ln2 = 5.72739- ■ ■ , p = pmax = 1- (b) The liquid-vapor 
coexistence curve of a conventional kagome lattice gas having nearest-neighbor pair interactions, 
with p being the particle number density and x = K^cjK^ being a reduced temperature where 
K^c = ln[(3 -|- 2\/3)] = 1.86626- - - . The coexistence curve is given by Pi°^^ = ^(1 ± m^), where 
nis is the spontaneous magnetization of a standard kagome Ising model. The rectilinear diameter 
of the symmetric rounded coexistence region is the horizontal solid line p = 1/2, which ends at 
the critical point (solid circle) whose coordinates are x = 1, p = pc = 1/2. The continuation 
of the rectilinear diameter into the disordered region is the horizontal solid line which continues 
indefinitely. 
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FIG. 5: (a) Elementary triplet correlation (012) vs reduced temperature x = K^c/K^{= T/Tc), 
where K^c = 3.96992 • • • . The solution (dashed curve) is determined along the coexistence curve 
including the critical point (sohd circle) with coordinates a; = 1, (012)c = \{27 + 15^3)7(26 + 
15-v/3) = 0.50961 • • • . The solution is also secured along the continuation of the curvilinear diameter 
(solid curve) into the disordered region, viz., the monotonically increasing and concave- upward 
curve, eventually ending at the point (solid square) with coordinates x = K3c/ln2 = 5.72739- • • , 
(012) = 1. (b) Conventional pair- interaction kagome lattice gas. Elementary triplet correlation 
(012) vs reduced temperature x = K2C/K2 where K2c = 1.86626 • • • . The solution (dashed curve) 
is determined along the coexistence curve including the critical point (solid circle) with coordinates 
X = 1, (012)c = ]^(3 + 2\/3) = 0.40400 • • • . The solution is also secured along the continuation 
of the curvilinear diameter (solid curve) into the disordered region, viz., the continuing sigmoidal 
curve which is monotonically decreasing and concave upward, eventually approaching (i)^ = 0.125 
at infinite temperature. 
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FIG. 6: Nine sites of the kagome lattice are specifically enumerated for selected applications of the 
theory. 
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FIG. 7: (a) Nearest neighbor pair correlation (01) vs reduced temperature x = K^c/K^{= T/Tc) 
where K^c = ln[(2 + \/3)^ + 1] = 3.96992 • • • . The solution (dashed curve) is determined 
along the coexistence curve including the critical point (solid circle) with coordinates x = 1, 
(01)c = 0.5239 • • • . The solution is also secured along the continuation of the curvilinear di- 
ameter (solid curve) into the disordered region, viz., the continuing sigmoidal curve which is 
monotonically increasing and slightly concave downward, eventually ending at the point (solid 
square) with coordinates x = ii'3c/ln2 = 5.72739- ••, (01) = 1. (b) Conventional pair- 
interaction kagome lattice gas. Nearest-neighbor pair correlation (01) vs reduced temperature 
X = K2c/K2{= T/Tc) where = ln[(3 + 2-v/3)] = 1.86626- - - . The solution (dashed curve) is 
determined along the coexistence curve including the critical point (solid circle) with coordinates 
X = 1, (01)c = n(i + \/3) = 0.43600- - - . The solution is also secured along the continuation 
of the curvilinear diameter (solid curve) into the disordered region, viz., the continuing sigmoidal 
curve which is monotonically decreasing and concave upward, eventually approaching (i)^ = 0.25 
at infinite temperature. 
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FIG. 8: Quartet correlation (1234) vs reduced temperature x = K3c/K3{= T/T^) where K^^ = 
3.96992 • • • . The sokition (dashed curve) is determined along the coexistence curve including the 
critical point (solid circle) with coordinates x = 1, (1234)c = 0.4617 • • • . The solution is also secured 
along the continuation of the curvilinear diameter (solid curve) into the disordered region, viz., a 
concave upward curve having an asymmetric rounded minimum before monotonically increasing 
towards an end point (solid square) with coordinates x = ii'sc/ In 2 = 5.72739 • • • , (1234) = 1. 




FIG. 9: Quintet correlation (01234) vs reduced temperature x = K^c/K-i{= T/Tc) where K^c = 
3.96992 • • • . The solution (dashed curve) is determined along the coexistence curve including 
the critical point (solid circle) with coordinates x = 1, (01234)c = 0.4596 ••• . The solution is 
also secured along the continuation of the curvilinear diameter (solid curve) into the disordered 
region, viz., a concave upward curve having an asymmetric rounded minimum before monotonically 
increasing towards an end point (solid square) with coordinates x = K-^d In 2 = 5.72739 • • • , 
(01234) = 1. 



32 



0.2- 1 

qI ^ i ^ 1 ^ 1 ^ 1 ^ 1 ^ 

"0 1 2 3 4 5 6 
X 



FIG. 10: Nonuplet correlation (012345678) vs reduced temperature x = K^c/K^{= T/Tc) where 
K^c = 3.96992 • • • . The solution (dashed curve) is determined along the coexistence curve including 
the critical point (solid circle) with coordinates x = 1, (012345678)c = 0.3970- • • . The solution 
is also secured along the continuation of the curvilinear diameter (solid curve) into the disordered 
region, viz., a concave- upward curve having an asymmetric rounded minimum before monotonically 
increasing towards an end point (solid square) with coordinates x = iC3c/ln2 = 5.72739- ••, 
(012345678) = 1. 
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FIG. 11: Probability p of an elementary cavity vs reduced temperature x = Ksc/K^{= T/Tc) where 
= 3.96992 • • • . The solution (dashed curve) is determined along the coexistence curve including 
the critical point (solid circle) with coordinates x = 1, Pc = (1234)c — (01234)c = 0.0021 • • • . The 
solution is also secured along the continuation of the curvilinear diameter (solid curve) into the 
disordered region, viz., a concave-downward curve with a symmetric rounded maximum located 
at X = 3.6 • • • , p = 0.0409 • • • , before monotonically decreasing towards a nodal end point (solid 
square) at x = K^d In 2 = 5.72739 • • • (fully occupied lattice). 
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FIG. 12: Enlargement of Fig 11 in the range of (reduced) condensation temperatures < a; < 1. 
Along the upper liquid branch pi of the coexistence curve, the probability p of an elementary cavity 
attains an asymmetric rounded maximum located at a; = 0.98 ■ ■ ■ , p = 0.0036 • • • . The solid circle 
is the critical point with coordinates x = 1, Pc = 0.0021 • • • . 
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